Search for effects beyond the Born approximation in polarization transfer observables 

in ep elastic scattering 
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Intensive theoretical and experimental efforts over the past decade have aimed at explaining the 
discrepancy between data for the proton electric to magnetic form factor ratio, Ge/Gm, obtained 
separately from cross section and polarization transfer measurements. One possible explanation for 
this difference is a two-photon-exchange (TPEX) contribution. In an effort to search for effects 
beyond the one-photon-exchange or Born approximation, we report measurements of polarization 
transfer observables in the elastic H(e, e!p) reaction for three different beam energies at a fixed 
squared momentum transfer — 2.5 GeV^, spanning a wide range of the virtual photon polariza- 
tion parameter, e. From these measured polarization observables, we have obtained separately the 
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ratio R, which equals UpGe/Gm in the Born approximation, and the longitudinal polarization trans- 
fer component Pe, with statistical and systematic uncertainties of AR ~ ±0.01(stat) ± 0.013(syst) 
and APi/Pf""'^ « ±0.006(stat) ± O.Ol(syst). The ratio R is found to be independent of e at the 
1.5% level, while the e dependence of Pe shows an enhancement of (2.3 ± 0.6)% relative to the Born 
approximation at large e. 



After decades of experimental and theoretical ef- 
forts, the internal structure of the nucleon remains one 
of the defining problems of nuclear physics. Based on the 
generally accepted notion that the electromagnetic inter- 
action is well understood from a theoretical point of view, 
elastic electron-nucleon scattering has served as a pow- 
erful tool to measure fundamental observables: the elec- 
tromagnetic form factors. There are two experimental 
methods for extracting the ratio of the electric to mag- 
netic form factors of the proton, Ge/Gm, from electron- 
proton elastic scattering. In the Rosenbluth separation 
technique [1 , G^ and G\j are determined from the an- 
gular dependence of the reduced cross section cr^ at con- 
stant Q^. Polarization experiments determine Ge/Gui 
by using a polarized electron beam with either a polar- 
ized proton target, or a measurement of the transferred 
polarization to the scattered proton. At values of the 
squared- momentum-transfer, < 1 GeV^, data from 
the Rosenbluth and polarization techniques are in good 
agreement. At large Q^, however, the cross section data 
PHS] disagree with the ratios Ge /Gm obtained using the 
polarization transfer method (SHH]- This systematic dif- 
ference is a source of intense debate in both the theoret- 
ical and experimental nuclear physics communities. 

Recently, theoretical attention has been paid to the set 
of radiative corrections (TUl |TT] that must be made to the 
cross section data in order to extract the form factor ra- 
tio. These corrections change the slope of the reduced 
cross section by as much as 30% for larger Q^. In con- 
trast, radiative corrections to the polarization data are 
essentially negligible [12] . Until recently, only "standard" 
radiative corrections were taken into account. Based on 
the observed discrepancy in the data, new efforts have 
been made to include higher-order radiative mechanisms, 
such as two-photon exchange (TPEX) [rSHTO] . Several of 
these calculations have indicated that TPEX partially re- 
solves the disagreement between the two data sets, but 
further investigation is needed. What is still lacking is 
a complete set of elastic ep scattering observables sensi- 
tive to the TPEX amplitudes, with sufficient precision to 
guide the development of a consistent theoretical frame- 
work for the interpretation of the discrepancy in terms of 
TPEX. This experiment is an effort to provide additional 
data. 

When considering the exchange of two or more pho- 
tons, the hadronic vertex function can generally be ex- 
pressed, in terms of three independent and complex am- 
plitudes, Ge,m = Ge,m{Q^) + SGE,M{Q^,e) and F3 
which are functions of and the kinematical parameter 
e = il + 2{l + r)tan2 0e/2)-\ where r = Q^/AAP, M 



is the proton mass, and 9e is the electron scattering an- 
gle. In the Born approximation, the first two amplitudes 
equal the real electric and magnetic Sachs form factors 
which depend only on Q^, and F3 vanishes. The exper- 
imental observables used to extract Gsp and Gmp from 
cross section and polarization transfer measurements, as- 
suming the validity of the Born approximation, are af- 
fected in different ways by TPEX, as shown in Eq.(l). 
The Rosenbluth method relies on measuring the e de- 
pendence of the reduced cross section cTj. at fixed to 
separate G"^ from G^j, and becomes highly sensitive to 
additive TPEX effects at large when tG\/TG\^ be- 
comes small. The transferred polarization to the recoil 
proton in H(e,e'p) has transverse (Pt) and longitudinal 
(Pi) components with respect to the momentum trans- 
fer in the scattering plane [501 HI] ■ The ratio R defined 
in Eq.(l) equals ^IpGe/Gm in the Born approximation, 
and is much less vulnerable to TPEX corrections. The 
TPEX corrections appear in ar, Pt, Pi, and R as interfer- 
ence terms between the Sachs form factors and the real 
part of the TPEX ampfitudes [H]: 
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where 3? stands for the real part, h — ±.1 and Pe are 
the helicity and polarization of the electron beam, and 

V — M ^"^^^ , with i?e, -Eg being the energy of the inci- 
dent and scattered electron, respectively. In the Born ap- 
proximation, these corrections vanish and the well known 
expressions for these observables [501 HI] are recovered. 
Other observables, such as the induced normal polariza- 
tion component, and the target-normal and beam-normal 
single-spin asymmetries, depend only on the imaginary 
(absorptive) part of the TPEX amplitude. A direct 
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test of the TPEX effect is the comparison between e~^p 
and e~p elastic scattering cross sections. Since the two- 
photon contributions (relative to the Born amplitudes) 
are of opposite sign, a few percent deviation from unity 
as a function of e is predicted for the ratio af,+ /cr^- . Re- 
cent analyses of e^p cross sections are inconclusive due 
to large uncertainties in the data [231^25] . 

In this experiment, carried out at Jefferson Lab in Hall 
C, a longitudinally polarized electron beam (82-86% po- 
larization) was scattered elastically off a 20 cm liquid hy- 
drogen target at = 2.5 GeV^. Electrons were detected 
by a 1744 channel lead-glass electromagnetic calorimeter 
(BigCal), which measured their coordinates and energy. 
Overlapping analog sums of up to 64 channels were used 
to form the BigCal trigger with a threshold of about half 
the elastic electron energy. Coincident protons were de- 
tected in the High Momentum Spectrometer (HMS) [25] ■ 
The HMS trigger was formed from a coincidence between 
a scintillator plane located behind the drift chambers and 
an additional paddle placed in front of the drift chambers. 
The polarization of scattered protons, after undergoing 
spin precession in the HMS magnets, was measured by 
the Focal Plane Polarimeter (FPP), which consists of an 
assembly of two 55 cm thick CH2 analyzer blocks, each 
followed by a pair of drift chambers to track re-scattered 
protons with an angular resolution of approximately 1 
mrad. Elastic event selection was performed offline in 
the same way as explained in [S], resulting in a very 
small inelastic contamination for all three kinematics; at 
e = 0.15, where it is the highest, the background fraction 
is 0.7%. 

The scattered proton polarization was obtained from 
the angular distribution of protons scattered in the ana- 
lyzer blocks of the FPP. The polar and azimuthal scat- 
tering angles {d, (p) of single-track events in the FPP 
chambers were calculated relative to the incident track 
defined by the HMS drift chambers. The difference 
and the sum of the azimuthal angular distributions for 
positive and negative beam helicities give the physical 
(helicity-dependent) and instrumental or false (helicity- 
independent) asymmetries at the focal plane, respec- 
tively. 

Since the proton polarization components are mea- 
sured at the focal plane, knowledge of the spin transport 
matrix of the HMS is needed to obtain and Pe at the 
target. The differential- algebra based modeling program 
COSY [27] was used to calculate the spin-transport ma- 
trix elements for each event from a detailed layout of 
the HMS magnetic elements. The quantities PeAyPt and 
PeAyPg were extracted using the maximum-likelihood 
method described in [S] [3] , with Ay the analyzing power 
of p + CH2 — ^ one charged particle -|- X scattering. Their 
ratio gives PtjPi independent of Ay and Pe- 

As an example of the quality of the data, Figjl] shows 
P and AyPg as a function of the vertical {dx/dz) and hor- 
izontal (dy/dz) slopes of the scattered proton trajectory 



relative to the HMS optical axis. Owing to the small 
acceptance of the HMS in both e and for all three 
kinematics, P and AyPi are constant across the accep- 
tance to a very good approximation. The absence of 
anomalous dependence of the extracted P and Ay Pi on 
the reconstructed kinematics is thus an important test of 
the accuracy of the field description in the COSY calcu- 
lations. In each panel of Fig. [l] the data are integrated 
over the full acceptance of all other variables. The hori- 
zontal line shows the one-parameter fit to the extracted 
data. In all panels, the per degree of freedom is close 
to one, indicating the excellent quality of the precession 
calculation. 
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FIG. 1: R (left column) and AyPt (right column) versus the 
dispersive dx/dz (vertical) and non-dispersive dy/dz (horizon- 
tal) slopes for the low-energy setting using the COSY model. 

The main results of this experiment are given in Ta- 
ble |T] and shown in Fig. [2] Figure [2^ displays P as a 
function of e with selected theoretical estimates. The 
data do not show any evidence of an epsilon dependence 
of P at = 2.5 GeV^. Both statistical and point-to- 
point systematic uncertainties (relative to the largest e 
kinematic) are shown in the figure. The total system- 
atic uncertainties in P are shown in Table |I] For a given 
data point, the point-to-point systematics are obtained 
as the quadrature sum over the differences between each 
of the systematic contribution and the corresponding one 
for a reference kinematic. Because the dominant sources 
of systematic uncertainty described below affect P for 
all three kinematics in a strongly correlated fashion, the 
systematic uncertainty on the relative variation of P as 
a function of e, characterized by the point-to-point un- 
certainties, is very small. 

Another sensitive probe of two-photon effects is taking 
the ratio of the measured Pi to Pf°^", where p^"^" is 
Pi calculated in the Born approximation. In the limit 
e — )■ 0, angular momentum conservation requires P^ — 1, 
independent of P (see Eq. [T]); for our measurement at 
e = .15, Pe varies by only 1.4% (relative) pS! for P be- 
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TABLE I: Kinematic table with the average quantities: the beam energy Ec , the momentum transfer squared , the electron 
scattering angle 9^, and the kinematical parameter e. Both the ratio R and longitudinal polarization Pt divided by Born 
approximation pp°^" are given with statistical (stat.), total systematic (tot.) and point-to-point (p.t.p.) uncertainties relative 
to the highest e point for R and to the lowest e kinematic for Pe/Pf°^". 
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tot. 
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tween and 1. Therefore, the measured value of hAyPi 
at e = .15 determines Ay = 0.15079 ± 0.00038 (specific 
to this polarimeter), corresponding to a relative uncer- 
tainty of 0.25%, included in the statistical error budget 
for Pf/Ff ''™. Applying the same phase space cuts at 
the focal plane results in Ay being the same for all three 
kinematics, at the lO^'^ level, p^"™ was calculated from 
the beam energy, the proton momentum, and the fitted 
value of R from this experiment, vifith the errors in each 
quantity accounted for in the total systematic error in 
Pi/Pf""'. In Fig. [2)3, the ratio Pg/Pf""^" is plotted ver- 
sus e. The results show an enhancement at large epsilon 
of 2.3 ± 0.6% relative to the Born approximation. 

The number of sources of systematic uncertainty is 
drastically reduced by the fact that the beam hclicity 
and the polarimeter analyzing power cancel out exactly 
in the ratio of polarization components. Consequently, 
the spin precession uncertainty is the dominant contri- 
bution. Since the central proton momentum was fixed 
across the three kinematics, the spin transport matrix is 
identical, resulting in small point-to-point systematic un- 
certainties. The error Ai/ifcend = ±0.5 mrad in the non- 
dispersive bend angle, due to uncertainty in the HMS 
quadrupole positions, represents the largest contribution. 
The uncertainty of ±0.1% in the absolute determination 
of the proton momentum, which has a negligible effect 
on the precession uncertainty in Pt/Pe, leads to a rel- 
ative uncertainty in the kinematic factor y^T{l + e)/2e 
that is roughly 1/3 of the relative uncertainty in Pt/Pg 
at the lowest e, but negligible at higher e. Errors in the 
dispersive bend angle, the beam energy, and the scatter- 
ing angle in the FPP give smaller contributions. The in- 
clusion of instrumental asymmetry terms obtained from 
Fourier analysis of the helicity-independent asymmetries 
in the likelihood function induces a negative correction to 
i? (|Ai?| < 0.013) and P^/P/^°™ (|APf/P/^°™| < 0.004) 
for each e value. A systematic uncertainty equal to half 
the false asymmetry correction was included in the to- 
tal systematic uncertainty on R and P^/P^"™. A 1% 
absolute systematic uncertainty (0.5% point to point) 
from the MoUer measurements of the beam polarization 
was added to the error budget of Pf/Pf In order to 
minimize systematic differences in the spin transport cal- 
culation among the three kinematics, cuts were applied 
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FIG. 2: a) _R as a function of e with statistical uncertain- 
ties, filled circles from this experiment and open triangle from 
[7]. The theoretical predictions are from: [14] (hadronic), [13] 
(GPD), [IS] (COZ and BLW) and [H] (SF) offset for clarity 
by -0.006 with respect to the fit. The one parameter fit result 
is: R = 0.6923 ±0.0058. b) Pi/PP°'''^ as a function of e. The 
point-to-point systematic uncertainties, shown with a band in 
both panels, are relative to the largest e kinematic in a) and 
relative to the smallest e kinematic in b). The star indicates 
the e value at which the analyzing power is determined. 



to the focal plane trajectories of the data for the two 
larger e points to match the smaller acceptance of the 
point at e = 0.15. The program MASCARAD [12] was 
used to compute "standard" radiative corrections to R. 
Small, positive corrections Ai? of 1.2x10"'^, 1.4x10"'* 
and 0.7x10^4 were found for e = 0.15, 0.63, and 0.78, 
respectively. The corrections to Pi / pp°™ were found to 
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be even smaller. The results shown in Table U do not 
include these corrections. 

The theoretical curves shown in Figure [2^ make widely 
varying predictions for the epsilon dependence of R. 
The hadronic model of Blunden et al. [H], where all 
the proton intermediate states are taken into account 
via a complete calculation of the loop integral using 4- 
point Passarino-Veltman functions |29j, shows a signifi- 
cant positive TPEX contribution at small e. The inclu- 
sion of higher resonances makes almost no difference |15j . 
On the contrary, the partonic model of Afanasev et al. 
[T5] , where the TPEX takes place in a hard scattering of 
the electron by quarks which are embedded in the nucleon 
through the GPDs, predicts a negative TPEX contribu- 
tion. A pQCD calculation of Kivel and Vanderhaeghen 
[TO] , which uses two different light front proton distri- 
bution amplitude parametrizations, one from Chernyak 
et al. (COZ) [30] and the other one from Braun et al. 
(BLW) Sr, presents a behaviour similar to the partonic 
model. The limit of applicability of the GPD and pQCD 
models is shown by the vertical dashed line on the fig- 
ure. The electron structure function (SF) based model 
developed by Bystritskiy et al. [T7], which takes into ac- 
count all high order radiative corrections in the leading 
logarithm approximation, does not predict any measur- 
able e dependence of R. The GPD, hadronic and pQCD 
models, while in good agreement with the available cross 
section data, predict a deviation of R at small e due to 
modification in Pt; not seen in the results presented here. 
Refering to Eq.Q, R is directly proportional to the Born 
value Ge/Gm, so all the theory predictions, which use 
a Ge/Gm value from [SHE], can be renormalized by an 
overall multiplicative factor. The enhancement seen in 
P(/Pf°™ is not predicted by any models. The behiavior 
of R at large e implies the same deviation of Pt from its 
Born value as the one observed in Pg/Pf°™. 

The high precision data presented in this letter add 
significant constraints on possible solutions of Eq. ([T]) for 
the real part of the TPEX amplitudes [32] . In this exper- 
iment, no e dependence was found in i?, suggesting that 
the 27 amplitudes are small or compensate each other in 
the ratio. The study of the non-linearity of the Rosen- 
bluth plot, a precise measurement of the single spin asym- 
metries and the determination of the j a^- ratio are 
essential to fully understand, quantify and characterize 
the two-photon-exchange mechanism in electron proton 
scattering. 
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